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^Sj \ Abstract 

We have carried out 2-D simulations of core-collapse supernova explosions. The 
frj \ local neutrino radiation field is assumed to have its maximum value either at the 

symmetry (polar) axis or on the equatorial plane. These lead to the prolate and 
oblate explosions, respectively. We find that the gain of the explosion energy in the 
q prolate explosion evolves more predominately than that in the oblate one when the 

total neutrino luminosity is given. Namely, the prolate explosion is more energetic 
^ \ than the oblate one. 

One of the authors (Shimizu et al. 2001) showed for the first time that globally 
anisotropic neutrino radiation produces more powerful explosion than the spherical 
neutrino radiation does. In our previous study (Madokoro, Shimizu, & Motizuki 
2003), we improved the numerical code of Shimizu et al. and demonstrated that the 
globally anisotropic neutrino radiation yields more energetic explosion than spatially- 
fluctuated neutrino radiation does. Together with the result of this paper, we conclude 
that the globally anisotropic (prolate) neutrino radiation is the most effective way of 
increasing the explosion energy among various types of explosions investigated in 
these studies. We discuss the reason for this. Our result is suggestive of the fact that 
the expanding materials of SN1987A is observed to have a prolate geometry. 
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1. Introduction 

Aspherical explosion is one of the key issues in the study of core-collapse supernovae. 
Observations have been suggested that the spherical symmetry is broken in several supernova 
explosions. For example, speckle and spectropolarimetry observations have shown that the 
spectra of core-collapse supernovae are significantly polarized at a level of 0.5 to 4% (Mendez 
et al. 1988; Jeffery 1991; Wang et al. 1996; Leonard et al. 2000; Wang et al. 2001; Hdflich 
et al. 2002). Recent Hubble Space Telescope images and spectroscopy have revealed that 
SN1987A has an axially symmetric geometry (Wang et al. 2002). Wang et al. (2002) also 
reported that the polarization represents a prolate geometry which has been fixed in time. 
They discussed that a jet-driven model based on asymmetry associated with neutrino flow 
(Shimizu, Yamada, & Sato 1994; Fryer & Heger 2000; Shimizu et al. 2001) may account for 
the observed asymmetries of SN1987A. Large values of observed pulsar kicks ( > 400 km s -1 ; 
see, e.g., Fryer, Burrows, & Benz 1998) are also a manifestation of asymmetry. Furthermore, 
it has been theoretically demonstrated that most simulations assuming spherical symmetry fail 
to yield robust explosions (e.g., Liebendorfer et al. 2001). These facts necessarily lead us to 
multi-dimensional simulations. 

Several groups have performed 2-D and 3-D simulations. It was shown (Miller, Wilson, 
& Mayle 1993; Herant et al. 1994; Burrows, Hayes, & Fryxell 1995; Janka & Miiller 1996; 
Mezzacappa et al. 1998; Fryer & Heger 2000; Fryer & Wallen 2002; Kifonidis et al. 2003) that 
convection in the neutrino-heated region behind the supernova shock or inside a nascent neutron 
star plays an important role to increase the explosion energy. 

Some groups also noticed the role of anisotropic neutrino emission from a protoneutron 
star. Janka and Monchmeyer (1989a,b) carried out 2-D simulations with rotation in order to 
explain the observed properties of SN1987A. When a protoneutron star rotates, centrifugal 
force will deform the core into an oblate form, and hence the neutrino flux distribution will 
preferentially be concentrated on the rotational (polar) axis. They concluded that the neutrino 
flux along the rotational axis could become up to three times larger than that on the equatorial 
plane. Later, Keil, Janka, & Miiller (1996) performed 2-D simulations without rotation but 
taking into account the convection inside the protoneutron star (Ledoux convection). They 
found that the neutrino flux on the surface of a protoneutron star shows a fluctuated anisotropy 
of ~ 3-4 %. Recently, Kotake, Yamada, & Sato (2003) performed simulations of the rotational 
core-collapse and obtained anisotropic neutrino heating, depending on their initial model with 
the degree of anisotropy in the neutrino heating rate from close to unity to more than a factor 
of 10. 

Shimizu et al. (1994, 2001) performed a series of systematic 2-D simulations in order 
to study the effects of anisotropic neutrino radiation on the explosion mechanism itself. They 
found that, when the total neutrino luminosity is given, only a few percent of enhancement in 
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the neutrino flux along the axis of symmetry is sufficient to increase the explosion energy by a 
large factor. They also found that this effect saturates around a certain degree of anisotropy. 
Furthermore, it was found that the effect of anisotropic neutrino radiation becomes more promi- 
nent when the total neutrino luminosity is low, while the difference between the anisotropic 
and spherical models is significantly diminished if the total neutrino luminosity is sufficiently 
high. Note that they proposed a rotation as one of the origins of anisotropic neutrino radia- 
tion; anisotropy in neutrino radiation may also be originated from convection inside the nascent 
neutron star, or asymmetric mass accretion onto the protoneutron star. 

Shimizu et al. (2001) considered only the case in which the maximum value in the 
neutrino flux distribution was located at the axis of symmetry. Janka & Keil (1997), however, 
pointed out that the neutrino flux can also be peaked on the equatorial plane. Following their 
previous study (Keil et al. 1996), they carried out 2-D simulations with both rotation and 
convection above and inside a nascent neutron star. It was found, at a very early phase of 
neutrino emission, that the neutrino flux was peaked at the rotational axis. However, they 
found that at a later stage of evolution Ledoux convection occurs strongly on the equatorial 
plane while suppressed near the rotational axis due to a combination effect of centrifugal force 
and convection. Consequently, the neutrino flux on the equatorial plane was enhanced over 
that along the rotational axis. 

In this paper, we are going to study the effects of the oblate neutrino heating, suggested 
by Janka & Keil (1997), on the explosion dynamics. We focus the shock position, the explosion 
energy, and the asymmetric explosion. Our purpose in this article, which follows the previous 
studies (Shimizu et al. 2001; Madokoro, Shimizu, & Motizuki 2003), is to investigate what kind 
of neutrino radiation is favorable for a successful explosion through a parameter study when 
the total neutrino luminosity is given. 

This paper is organized as follows. In § 2, we describe our numerical simulation. We 
then discuss our numerical results in § 3. The results of the oblate explosion are compared 
with those of the prolate one. We shall find that the locally intense neutrino heating along the 
pole, which leads to a prolate explosion, is more effective way of increasing the explosion energy 
than the oblate heating when the total neutrino luminosity is given. Finally, a brief summary 
is presented in § 4. 

2. Numerical Simulation 

We perform 2-D simulations by solving hydrodynamical equations in the spherical co- 
ordinate (r,8). A generalized Roe's method is used to solve the general equations of motion. 
Our computational region ranges from 50 km to 10000 km in radius from the center of the 
protoneutron star. We start our calculations when a stalled shock wave is formed at ~ 200 km. 
The temperature on the neutrinosphere, T u , and the total neutrino luminosity, L u , are assumed 
to be constant in our simulations. The details of our numerical technique, the EOS and the 
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initial condition, are described in Shimizu et al. (2001). 

We obtain the initially stalled shock by solving stationary hydrodynamic equations with 
assuming spherical symmetry For the purpose of this paper, this assumption is enough. Note 
that in a simulation model by other group (Kotake et al. 2003), the stalled shock wave can 
become almost spherical even if the protoneutron star rotates with the rotational velocity which 
the authors call 'moderate' (and with weak differential rotation). We stop our calculations with 
constant L v at t = 500 ms after the shock stagnation. This is because recent one- and two- 
dimensional simulations have shown that the timescale of decay of the neutrino luminosity after 
the shock stall is about 500 ms (Bruenn 1993; Woosley et al. 1994; Janka & Miiller 1995; Fryer 
& Heger 2000). 

We fix the radius of the protoneutron star, -Rns ; to be 50 km in our simulation. The 
neutrinosphere actually shrinks as a result of neutrino emission (Burrows et al. 1995). However, 
the effect of shrinking of the neutrinosphere is not important in our present study. The reasons 
are as follows: Firstly, the degree of anisotropy in the neutrino radiation at a distant point from 
the neutrinosphere (which we denote 'C2') depends strongly on the degree of anisotropy on the 
neutrinosphere (denoted by 'a'), while is insensitive to -Rns (see equation (8) later). Secondly, 
the region in which we are interested is 1000 km above the neutrinosphere. We then expect 
that small change of the radius of the neutrinosphere (from 50 km to 40 km; see Burrows et 
al. 1995) does not affect the shock dynamics qualitatively during the timescale of 500 ms on 
which we simulate. Thirdly, we fix both the radius of the protoneutron star and the neutrino 
energy (or equivalently the neutrino temperature) in our simulation. On the other hand, the 
neutrinosphere actually shrinks as a result of neutrino emission, while at the same time the 
neutrino energy increases. It is therefore a good approximation to consider that the rate of 
absorption of the neutrino energy on the shocked matter is approximately constant during the 
short timescale of 500 ms. We then think that the qualitative result is the same as our present 
result, even if we include both the shrinking of the protoneutron star and the increase of the 
neutrino energy. 

We have improved (Madokoro et al. 2003) the numerical code of Shimizu et al. (2001): 
the cells in the ^-direction were shifted by half of the cell size (Shimizu 1995) in order to avoid 
a numerical error near the pole. Note that this numerical error was not serious and minor for 
the investigation of the explosion energy, but may affect the results of nucleosynthesis. 

The local neutrino flux is assumed as 



where a is the Stefan-Boltzmann constant. In equation (1), c 2 is a parameter which is related 
to the degree of anisotropy in the neutrino radiation. In order to see the effect of anisotropic 
neutrino radiation itself on the explosion, the value of C\ is calculated from given c 2 so as to 
adjust the total neutrino luminosity to that in the spherical model at the same T v . The total 
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neutrino luminosity is obtained by integrating equation (1) over the whole solid angle, 

L V = J r%(r,9)dQ = ^aT* 4tt Ci (l + ^c 2 ) , (2) 



which is equated to that of spherical explosion with the same Tj 



L? = Y^Tt 4„R* S . (3 ) 
By comparing equations (2) with (3), we obtain 

One can easily confirm that the neutrino flux along the pole (Z 2 ) and that on the equa- 
torial plane (l x ) are proportional to ci(l + c 2 ) and ci, respectively. The degree of anisotropy at 
a distant point far from the neutrinosphere, l z /l x , is then given by 

l f = l + c 2 . (5) 

''X 

Accordingly, the neutrino flux has its maximum at the pole or on the equator when the sign of 
c 2 is positive or negative, respectively. 

We have to emphasize that the degree of anisotropy in the neutrino flux distribution 
for an observer far from the neutrinosphere is different from that on the neutrino-emitting 
surface. This is schematically illustrated in Fig. 1. When we observe the neutrino flux far from 
the neutrinosphere, the local neutrino flux is seen as equation (1). On the other hand, the 
neutrino flux on the neutrino-emitting surface has a similar profile but the degree of anisotropy 
is different from equation (1). This is because a distant observer obtains the neutrino flux by 
integrating all the contributions over the solid angle from an anisotropically radiating surface, 
and hence the degree of anisotropy is reduced for the observer (the geometric effect we call). 
For the neutrino flux on the neutrino-emitting surface, c 2 in equation (1) is replaced by a, 
where a is a parameter which represents the degree of anisotropy in the neutrino radiation on 
the neutrinosphere. The angular dependence of the neutrino flux on the neutrinosphere is then 
represented as 

l v ~ l + acos 2 #. (6) 

In principle, the value of c 2 is calculated from a given a taking into account the geo- 
metric effect from an anisotropically radiating surface. Although it is difficult to calculate the 
exact relationship between c 2 and a, we can estimate it in a same way as in Madokoro et al. 
(2003). The neutrino flux observed far from the neutrinosphere is obtained by averaging all 
contributions from the flux on the surface of the neutrinosphere. From the similar procedure 
to that in Madokoro et al. (2003), we obtain the ratio of the local neutrino flux along the polar 
axis (l z ) to that on the equatorial plane (l x ) for an observer far from the neutrinosphere as a 
function of a, 
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^i±^ ( 7) 

l x 4 + a ' K J 

By comparing equation (7) with equation (5), we finally obtain 

c 2 ~— . (8) 

Note that equation (7) is different from equation (7) in Madokoro et al. (2003). This is because 
the profile of the neutrino flux on the neutrinosphere is approximated by a step function in 
Madokoro et al. (2003), while we directly use equation (6) to obtain (7) in this paper. 

3. Prolate and Oblate Explosions 

We first study the prolate and oblate explosions with the same degree of anisotropy; 
l z /l x = 1.1 (model prolOO) and l z /l x = 1/1.1 (model obl091). The value of c 2 is 0.1 and -0.091, 
respectively. The spherical model (model sphOOO) is also examined for comparison. The temper- 
ature on the neutrino-emitting surface, T v , is at first fixed to be 4.70 MeV. For models prolOO, 
obl091, and sphOOO, we also study the cases with different values of the neutrino temperature 
(T v = 4.65 and 4.75 MeV) in order to check the sensitivity of the results on the neutrino tem- 
perature. We further consider several cases changing the degree of anisotropy in the neutrino 
radiation in order to examine at what degree of anisotropy the effect of anisotropic neutrino 
radiation will saturates for both the prolate and oblate models. The models we examined are 
summarized in Table 1. 

Figure 2 depicts the color-scale maps of the dimensionless entropy distribution with the 
velocity fields for the two models, prolOO and obl091, at t ~ 280 ms after the shock stall (for 
translation to the dimensional entropy, see Shimizu et al. 2001). The color boundary between 
dark blue and light blue shows the shock front at r ~ 2500 — 3800 km for model prolOO and 
at r ~ 1800 — 2200 km for model obl091. In model prolOO, higher entropy upflows are mainly 
formed in the region of high latitudes due to intense neutrino heating on the polar axis. As the 
locally heated matter pushes the shock front, the shock becomes prolate and largely distorted. 
On the other hand, the entropy distribution of model obl091 shows that the shock front is 
elongated in a oblate shape as a result of neutrino heating concentrated on the equatorial 
plane. We also notice that the shock front in model prolOO is clearly more extended than that 
in model obl091. 

The evolution of the explosion energy is illustrated in Fig. 3. The result of the spherical 
model is also shown for comparison. We find that the energy gains are substantially different 
between the prolate and oblate models. The increase in the explosion energy is much larger in 
model prolOO than that in model obl091, although both models finally appear to explode. 

The difference between the prolate and oblate models is attributed to the mechanism 
which was firstly pointed out by Shimizu et al. (2001): importance of the locally intense neutrino 
heating. It should be noted here that the situations of prolate and oblate explosions are different. 
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In the former, neutrinos are predominantly emitted along the polar axis. Because of increased 
pressure in the locally heated matter near the pole, the shock wave is partly pushed outwards in 
the polar direction. The other part of the shock will follow due to the pressure gradient along 
the shock front (Rankine-Hugoniot relation). In Fig. 3, one sees that the explosion energy 
begins to increase at the position of t ~ 40 ms. It is confirmed that this point corresponds to 
the time when the stalled shock starts to move outwards. 

On the other hand, in the case of the oblate explosion, the neutrino radiation is 
anisotropic but its peak is not located point-like; the heated region is disk-like over the equa- 
torial plane. The efficiency of anisotropic neutrino heating of the oblate model is considerably 
reduced compared with that of the prolate (local) one. This is naturally expected from the 
physical mechanism claimed in Shimizu et al. (2001). 

Note that the neutrino heating rate itself is proportional to and there are no re- 
markable difference in the neutrino heating rate between the prolate and oblate models. The 
neutrino cooling rate, however, behaves as ~ T^, where T m is the matter temperature. As a 
result of the earlier shock revival, T m in model prolOO is quickly reduced. The neutrino cooling 
for the prolate model is therefore largely suppressed as the decrease of matter temperature in 
the course of the shock expansion. Thus, the explosion energy in model prolOO becomes larger 
than that in model obl091 (see Fig. 3). 

Up to here, we have shown only the results of two models (prolOO and obl091) which 
have the same degree of anisotropy in the neutrino radiation. We need not to study many 
models for the prolate and oblate explosions with the same degree of anisotropy. From sys- 
tematic investigation in the previous papers (Shimizu et al. 2001; Madokoro et al. 2003), it is 
learned that we will obtain the same result if we examine many models with various degrees of 
anisotropy for the prolate and oblate neutrino radiation: the prolate explosion is always more 
powerful than the oblate one with the same degree of anisotropy for a given luminosity. 

Next, we investigate the cases in which the prolate and oblate explosions have different 
degree of anisotropy. We have found that the effect of anisotropy becomes more powerful in both 
the prolate and oblate models as the degree of anisotropy in the neutrino radiation increases: 
the larger the degree of anisotropy is, the larger the explosion energy becomes. It was found, 
however, that this effect saturates around a certain degree of anisotropy (see discussion in 
Shimizu et al. 2001 and Madokoro et al. 2003). Figures 4 and 5 show the saturation properties 
of the anisotropic neutrino radiation for the prolate and oblate models, respectively. We notice 
that, at a later stage of evolution, model prolOO (t < 360 ms) or prol50 (t > 360 ms) shows the 
largest explosion energy among the prolate models (Fig. 4), while the explosion energy of model 
obll67 becomes largest among the oblate models (Fig. 5). This means that the saturation of 
the effect of anisotropic neutrino radiation comes at l z /l x ~ 1.1 — 1.15 for the prolate explosion 
and l x /l z ~ 1.2 for the oblate explosion. We also confirmed that the explosion energy in the 
prolate model of l z /l x = 1.1 is larger than that in the saturated oblate model of l x /l z = 1.2. This 
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is explicitly illustrated in Fig. 6. In other words, the prolate explosion of l z /l x — 1.1 is always 
more energetic than any other oblate models even if the degree of anisotropy in the neutrino 
radiation becomes larger. We therefore conclude that the prolate explosion is generally more 
energetic than the oblate one for a given neutrino luminosity. 

For the purpose of clarifying the sensitivity of the supernova problem on the neutrino 
heating, we have also carried out simulations of less luminous [T v = 4.65 MeV) and more 
energetic (T u = 4.75 MeV) models. When T v = 4.65 MeV, we found that the prolate model 
(prol00T465) does explode while the oblate one (obl091T465) fails. Note that the difference 
between 4.70 and 4.65 MeV is only 1%. In contrast, the difference between the prolate and 
oblate models becomes much smaller when T u is increased to 4.75 MeV (models prol00T475 
and obl091T475). These results indicate that the effect of local neutrino heating becomes more 
pronounced and therefore more important when the total neutrino luminosity is lower. 

The spherical model also finally explodes when T v = 4.70 MeV or higher (see Fig. 3). 
This simply means that any models will explode irrespective of the degree of anisotropy if 
the luminosity is sufficiently high. However, we should keep in mind that the total neutrino 
luminosity can not be simply increased due to the problem of small mass of the protoneutron star 
and that of Ni overproduction, especially in the case of essentially spherical models (Burrows 
et al. 1995; Janka et al. 1996). On the other hand, our model based on anisotropic neutrino 
radiation can explode at lower neutrino luminosity, and therefore could solve the problems 
described above. 

In our previous article (Madokoro et al. 2003), we made a comparison between globally 
anisotropic (prolate) explosions and fluctuated explosions. We concluded in that paper that 
globally anisotropic neutrino radiation is more effective than the fluctuated neutrino radiation 
to increase the explosion energy for a given neutrino luminosity. In addition, in the present 
study we have found that the prolate neutrino heating produces more energetic explosion than 
the oblate one. Combined with the result of Shimizu et al. (2001), we conclude that the globally 
anisotropic (prolate) neutrino radiation is the most effective way of increasing the explosion 
energy among various types of explosions investigated in these studies when the total neutrino 
luminosity is given. These results support the statement made by Shimizu et al. (2001): globally 
anisotropic neutrino radiation, or locally intense neutrino heating is of great importance to 
produce a successful explosion. 

As we have shown, the locally intense neutrino heating along the polar axis produces a 
robust explosion in which the shock front is deformed in a prolate form. Our results may be 
related to the observation that the expanding materials of SN1987A has a prolate geometry 
(see Wang et al. 2002, Hoflich et al. 2002). 
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4. Conclusion 



We have performed 2-D numerical simulations of core-collapse supernova explosions with 
prolate and oblate neutrino radiation fields. We found that these two models give different 
results in the shock position at the same evolutionary stage of explosion, the explosion energy, 
and the geometry of asymmetric explosion. It is also found that the prolate explosion yields 
larger energy gain than the oblate one for a given neutrino luminosity. We can expect from the 
mechanism suggested by Shimizu et al. (2001) that the prolate explosion is more energetic than 
the oblate one when the degree of anisotropy in the neutrino radiation is the same between 
these two. It is not obvious, however, whether the prolate explosion is always more powerful 
than the oblate one when the degree of anisotropy varies independently in these two models. 
Our numerical study in this paper has clarified that the prolate explosion (l z /l x = 1.1) is always 
more powerful than the oblate ones irrespective of the degree of anisotropy. 

Moreover, we found that the difference between the prolate and oblate explosions be- 
comes prominent when the total neutrino luminosity is low. This means that the local (prolate) 
neutrino radiation becomes more important to increase the explosion energy as the total neu- 
trino luminosity decreases. Because we can not simply increase the total luminosity to explain 
the observed explosion energy due to the problem of small mass of the protoneutron star and 
that of Ni overproduction, especially in the case of essentially spherical models (Burrows et al. 
1995; Janka et al. 1996). Therefore, we conclude that the local (prolate) neutrino radiation is 
of great importance in actual supernova explosions. 

One of the authors showed for the first time that the locally intense neutrino heating 
along the axis of symmetry produces more powerful explosion than the spherical explosion when 
the total neutrino luminosity is given (Shimizu et al. 2001). Our results in this article, together 
with the results of our previous work (Madokoro et al. 2003) show that the prolate explosion 
is the most effective way of increasing the explosion energy when the total neutrino luminosity 
is low among various types of explosions examined in these studies. Our results are suggestive 
of the observation that the expanding materials of SN1987A has been observed to be deformed 
in a prolate form (see Wang et al. 2002, Hoflich et al. 2002). 
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Fig. 1. A schematic picture which shows the relationship between the neutrino flux on the neutrinosphere 
and that for an observer far from the neutrinosphere. Due to the geometric effect from an anisotropically 
radiating surface, | c 2 1 is always smaller than \a\. 
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Fig. 2. Color-scale maps of the dimensionless entropy distribution and the velocity fields for the two 
models, prolOO (Left) and obl091 (Right) at t ~ 280 ms after the shock stall. 
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Fig. 3. Evolution of the explosion energy (E cxp \) for model prolOO (solid line) and model obl091 (dashed 
line). The result of the spherical explosion, model sphOOO (dotted line), for the same T„ is also shown for 
comparison. 
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Fig. 4. Evolution of the explosion energy (E exp \) for the prolate models: pro050 (dotted line), prolOO 
(solid line), prol50 (short-dashed line), pro200 (long-dashed line), and pro250 (dot-dashed line). The 
result of model pro300 is close to that of model pro250 and hence we omit it. 
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Fig. 5. Evolution of the explosion energy (-Eexpi) for the oblate models: obl048 (dotted line), obl091 
(short-dashed line), obll30 (long-dashed line), obll67 (solid line), and obl200 (dot-dashed line). The 
result of model obl231 is close to that of model obl200 and hence we omit it. 
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Fig. 6. Evolution of the explosion energy (-E e xpi) for the saturated models: prolOO (solid line), 
pro!50(dot-dashcd line), and obll67 (dashed line). 
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